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Q>^ • Under the name prime decomposition (pd), a unique decomposition of an arbitrary A^- 
^ , dimensional density matrix p into a sum of seperable density matrices with dimensions 
Q-i! determined by the coprime factors of A^ is introduced. For a class of density matrices 
^ _ a complete tensor product factorization is achieved. The construction is based on the 

■ Chinese Remainder Theorem, and the projective unitary representation of Zj^i by the 
^ ■ discrete Heisenberg group H^. The pd isomorphism is unitarily implemented and it is 

shown to be coassociative and to act on as comultiplication. Density matrices with 
complete pd are interpreted as grouplike elements of H^. To quantify the distance of 
^ I p from its pd a trace-norm correlation index £ is introduced and its invariance groups 
c3 ■ are determined. 
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Quantum correlations, an emblematic notion of quantum theory, remains an open 
challenge since the early days of Quantum Mechanics [^, Recent investigations have 
set important questions concerning classification of various types of quantum correla- 
tions and their appropriate quantification. These theoretical activities have parallel 
developments with, and are partly motivated from recent interesting proposals which 
engage quantum correlations to such diverge tasks as e.g quantum computation and 
communication quantum cryptography 0, teleportation 0, and some new fre- 

quency standards 0. Although the classification of quantum correlations is still open 
to refinements, it appears to include the following cases: for pure states, correlations 
entail nonlocality and give rise to violation of Bell inequalities [§] . For mixed states, two 
systems are considered uncorrelated if the composite system density matrix factorizes 
into a product of reduced density matrices, one for each isolated quantum subsystem 
viz. p = pi® p2, where pi^2 = Tri 2P, are determined by partial tracing. Quantifica- 
tion measures for that case include the von Neumann entropy and other invariant 
indices 0. On the other hand classical correlations for quantum subsystems imply 
seperability of the joint system density matrix, which is analysed into a convex sum 



for products of pure states viz. p = J^iPiPi ® P2-, ^ Pi < 1; Hi Pi = 1; [113 • Neces- 
sary and sufficient conditions for the existence of such convex decompositions for p's 
acting on C x C and C x C , became available recently 111], |12[. Upper bounds for 
the number of terms in such convex expansions of seperable matrices have also been 
determined, together with construction algorithms for the cases dimTi < 6 IT^ and 



dimH < oo WM. Beyond these types of classical correlations we encounter inseperable 



or entangle quantum states. For their characterization and the quantification of their 
entanglement some general conditions have been presented that good entanglement 
measures should satisfy |I5 |. 

In this Letter we address the problem of the prime decomposition (pd), of a finite 
but otherwise arbitrary A^-dimensional square density matrix p, into a sum of products 
of elementary density matrices, the number and the respective dimensions of which are 
determined by the compositeness of the dimension of p. This is achieved by means of 
1) the so called Chinese Remainder Theorem (CRT) |T^, that is based on the prime 
decomposition of (this also explains the name we have chosen for the decomposition), 
and 2) by the fact that the discrete Heisenberg group Hj^, provides a projective repre- 
sentation of the abelian cyclic group Zjv More concretely, if = p^^p^^ . . ■pT'', 
is the prime factor decomposition of N, where jo's are distinct primes, then the pd of 
the density matrix involves square matrices p^*-* i = 1, . . . , t, with power prime dimen- 
sion equal to Ni = Also the number t of p- factors is bounded by the number of 
coprime factors of N. As a measure of the correlation of a given mixed state p, with its 
possible prime or other decomposition, we evaluate the trace-norm distance between 
the two densities, study its unitary invariant symmetries, and interpret it in terms of 
the quantum variances between local operators of the subsystems of the decomposition. 

We start by considering the matrix realization of the discrete Heisenberg group 
generated by the operator set of A^ elements Jm = Jmim2 = ui2"^^'^^g'^'^h"^^, where the 
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matrices 



g = dmg{l,uj,...,u 
h = \n){n+l\ , (1) 



nez 



AT 



satisfy the relations ugh = hg, h'^ = h ^ , g^ = h'^ = 1, hh'^ = Wh = 1, gg'' = g'^ g = 1, 
with a;^ = 1, and (mi, 7712) G Z^, the square index-lattice. By virtue of these relations 



the following commutators are valid |]T8[, 
[Jm, Jn] = -2zsin 



TV 

— m X n 

N 



Jm+n (mod A^) . (2) 



Moreover due to linear independence, completeness and the orthonormality issued by 
the relation 

TrJmJn = N6m+n,o (mod N) , (3) 



the same generator set forms a basis of the su{N) matrices JTT 



Let us consider a A^-dimensional quantum system S with Hilbert space Hn- The 



generators of the finite Heisenberg group provide an operator basis {Jm\m G 
for the decomposition of the density matrix p of S, viz. 



-AT J 



P=^ ('^m^m) = ^[l+ ^n,Jm] ■ (4) 

ttigZjy meZjv 

with Z^ = Zjv X Zat \ (0,0), We note here that due to the Hermitian conjuga- 
tion of the basis elements i.e J}^ = J-mi-m2 = JN-mi,N-m2i the Hermiticity of 
the density matrix imphes for its elements the reality conditions \*^ = XN-m- Let 
us assume to be a composite positive integer with prime-power decomposition 

= Pi^P2^ ■ ■ -Pt* = N1N2 . . . Nt, where each of the factors is distinct, uniquely de- 
termined and relatively prime to each other, i.e gcd{Ni, Nj) = 1 when i ^ j. Then 
according to the CRT the isomorphism Z^r = Z^Vi © ■ ■ ■ © Z^Vt, is valid for the index- 
cyclic groups labelling the operator basis. To proceed we introduce the group iso- 
morphic map Z^ Z^^ © ■ ■ ■ © Z^^ , between the cyclic groups. The explicit defi- 
nition reads: (mi, 7722) (5(mi), 5(m2)) (mn; m2i, mi2; m22, • • ■ , m^; m2i), where 
rriii = nil — pNi, m2i = m2 — qNi, z = 1, . . . , t, p, g G Z, stand for the residues of the 
division of mi, m2 by Ni. 

Next we regard the fact that H^, provides a projective unitary representation 
of the additive cyclic group Z^, by means of the map Z^ -5 H]\r. More explic- 
itly, (mi,m2) ^ 7rAr(mi,m2) = Jmi,m2, with the property 7rAr(m + n) = J^n+n = 
JmJn^^"^^"' = 7r7v(m)7r7v(n)e2'^^", where m x n := min2 — m2ni. Then the following 
commuting diagram: 

Z^ ^ ^Ni © ■ • • © 

TTat i i TTjy^ X ■ ■ ■ X TTNt 

Hn — > ifjvi ® ■ ■ • © Hn, 

T^5 
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Fig. 1 

given by the equation ttsottn = (vtat^ x • • • x vTivJ o S, induces the isomorphism of CRT 
from the index-groups to the associated Heisenberg groups by the following component 
version of the above diagram: 

m S{m) ^ {mn; 17121,..., mit;m2t) 

TTAT i i TTjVi X • • • X TTjVt 

Fig. 2 

We state the main proposition for the prime decomposition: 

Proposition. The isomorphism tt^, determined by the commuting diagram of Fig. 
1, via its component version Fig. 2, is a linear map which induces the 5-map of CRT 
into the Heisenberg group i/jv, and provides the unique pd of elements of H^. Also 
TTs is implemented by unitary operator in the Hilbert space 71^ and possesses the 
coassociativity property. 

Proof. If m e and Jm = a;i/2'»i™2^mi/jm.2^ then 5{mi,m2) = {mu;m2i, . . . , mu, m2t), 
with mii and m2i, the residues of the division of mi,m2 by Ni respectively. Accord- 
ing to CRT 5 is an isomorphism the determination of which provides the solution 
of a system of congruences mi = mu (mod Ni) and m2 = m2i (mod iVj), when 
gcd(A^j, Nj) = 1, Ni ^ Nj, i.e when N^, Nj i = 1, . . . ,t, are pairwise coprime posi- 
tive integers. Inversely, given the residues, the numbers mi,m2 can be determined in 
a mixed-radix notation by mi = Yl\=i''^iiN iUi and m2 = J2l=i'>^2iNiyi, (mod N), 
where :— ^ and i/i is the solution of the congruence N^yi = 1, (mod Ni). Alter- 
natively by means of the Euler function (t){k), which counts the positive integers / < k, 

which are coprime to k, the yi is given by yi = N^^^^^ ^, (mod N^). Then mi,m2, 

are expressed in the form mi = J2l=i i^T'iiNf^^'^ and m2 = I]*=i m2iNf^^'^ (mod N). 

We turn now to study the consequences of this decomposition for the generators of 
Hisf. With the notation as before we obtain the relations 



, =n^r (5) 

i=l 



where Qi :— g^i't'(^i) and gj^^ — g^i^i<t>{^i) — g^H^t) — i. Analogous relations hold for 
the generator h"^'^. By direct computations it is verified that gihj — hjgi ii i ^ j, and 
gih\ — ujfh\gi, where Ui :— u^* ' . This definition implies that uji is periodic with 

respect to the coprime factors of N i.e LOi — a; * * — 1, for i — 1, . . . ,t. Using 
the above commutation properties of the generators we write: 



7r5 7r^(mi,m2) = 7r,( J^) = n ^^1^'™^^^/^^ ^ 11 •^S.m. 

1=1 



^i=iJmum2i = (^^1 X • • • TTjvJ (miim2i, . . . , mum2t) ■ (6) 
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The isomorphism introduced above is based on the fact the the Jm^^rn2i^ commut- 
ing for different i's and their moduh make them to behave as copies (vr^-isomorphic 
images ) of the original Jmim2 ^ Hjy, with periodicities Ni < N; this is similar to 
harmonics in Fourier analysis. The following embedding provides the explicit form of 
the isomorphism: 

with mii,m2i € Z^.; this provides the commutativity of the diagrams. 

The prime decomposition of a general density matrix given in eq. with coeffi- 
cients pm = Tr Ar(Jmp), can now be evaluated and reads, 

^s{p) ~ T7 5^ P5{mi,m2)^s{'Jmi,m2) — "77 P5{mi,m2)J&{mi,m2) 

(mi,m2)e^jv (mi,m2)e^jv 

= T7 5^ ■ ■ ■ Pmiim2i,...,mitm2t'^miim2i ® ' ' ' ® Jmum2t 

--w2 -y2 

(miim2i)e^iVi (mitm2t)G^jv4 

where Pmiim2i,...,mitm2t ~ Tr^Vi ■ ■ ■ Tr7Vt('7r5(/9) Jmiim2i ® ■ ■ ■ ® Jmitm2t)- 

This suggests that we first map p into T^sip), according to the previous analysis and 
then project along the J^^. ^^.'s , in order to determine the coefficients of the p-matrix 
factors in the prime decomposition. 

A special form of pd that contains only a single product term is possible for a 
special class of density matrices with coefficients pm = where /(m) G /2(Z^), 

an arbitrary real function. If /(m) = X]fe« /a:/'^i"12, the tranformed coefficients pum) = 
due to io^ = 1, factorize as follows: 



P5{m) 



N 
N 



= Il^,U^I''^''''^''^=--tlPin.u,rn2.). (9) 
i=l i=l 

Note that power raising is counted (mod A^), so above a Frobenious type of map 
has been used i.e m\ = (I]i=i 'muNf^^''^)^ = mjjivf*'^*'' (mod N). Therefore we 
obtain the pd ns{p) = p^^^ ® ■ ■ ■ ® p^*\ where p^^^ = P{mu,i^2i)J^iira2r In 

view of the coassociative property of the tt^, to be established shortly, we see that those 
matrices that admit such complete factorization behave as grouplike elements under 
the pd map. 

To proceed we study the unitary implementation of the pd of density matrices. We 
introduce the operator Vs : TIn — ^ ®\=iH-Ni, given by 

Vs= E IK})(^I= E \ni)®---®\nt){n\ , (10) 

nsi^N ngZjv n(i^isi 
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and its conjugate 

{ni}G{Z]v.} {njG{Zjv.} 

(11) 

(this S as the one used earher maps numbers to their respective residues; see below). 
These operators form a conjugate pair that obeys the unitarity condition VsVg = 
®i=ilwjv. '^s'^s = l-Hjv- Then it is straightforward to verify that the pd map ns act- 
ing on general density matrix is implemented by the unitary similarity transformation 

Finally we study briefly an important property of the pd map tt^, namely that it 
becomes a coassociative comultiplication of the Heisenberg group H]^\ this illustrates 
a connection of CRT with Hopf algebras in the framework of quantum mechanical 
correlations. Consider the map J„j^„2(^) = {x — pni,x — (Jn2), p,a & 'Z, by which a 
X G Z„^„2, decomposes into its resedues wrt coprimes ni,n2. Also consider its dual 
map p.ni,n2{ci-,b) = an2^^^^ + bnf^^^^ = x (mod nin2), which constructs the solution 
of the congruences x = a (mod rii), x = b (mod 71,2), according to CRT. Then we 
check that for Ni, N2, three coprime factors of A^, the following equation is valid on 
any a G Zat: 

i^NuNi X id) o 6niN2,N3 = {id x 5n2,N3) ° Sni,N2N3 ■ (12) 
This is dual to the relation 

PniN2,Ns ° {pni,N2 X id) = Pni,N2N3 ° {id x Pn2,N3) , (13) 

which holds if we are given three congruences and combine them pairwise in two dif- 
ferent ways. This (co) associativity of the CRT maps, in turn is induced into the pd 
map TTs, where it takes the form 

{^Sj,^,N2 ® id) o '^s^,N2.N3 = {id ® 71-5^2. ) ° ^-S^i.iv^iva " (14) 
As an example we take the system 

X = 2 (mod 3) , 
X = 2 (mod 4) , 

X = 3 (mod 5) , (15) 

with solution x = 38 (mod 60), and obtain 

/i3-4,5 o {P3A X id){2, 2, 3) = /i3.4,5(2, 3) = 38 

/i3,4.5 0(^(^x^4,5) (2, 2, 3) = /i3,4.5(2,18) = 38. (16) 

Dualizing we recover the relation for the 5's which induces the coassociativity of the 
pd mapping: 

{id O Tis^J o 7rs^^^,,{p) = {ivs^^^ ® id) o irs^,^^,{p) , (17) 
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for p G Hqq. Closing this proof we note that the integral • — ^ ^ with defi- 
nition Jj^ p := TtnP, is invariant under the comultiplication nsj^_^ in the sense that 

We turn now to the study of the correlation between finite quantum systems. We 
start with two systems with state vector Hilbert spaces of dimension Ni,N2 respec- 
tively. Any observable and density matrix is expressed by the elements of the Lie 
algebra u{Ni), u{N2) correspondingly. For the density matrix of system-1 e.g , 



.(1) 



iVi 



and similarly for system-2. The choice of the operator basis (l*-*-*, 



1,2), for 



the Lie algebra u(Ni) u{l 
the density matrix reads 



su{Ni) is an important one. For a composite system 



P 



+ 



E 



1(1)^1(2)+ y: A«J«®1(2)+ J2 Afl«®jf 

V A(1'2) 7(1) (9) 7(2)" 



-N2 



(19) 



where A« ^ (j^i)) = Tr(p • jW ® l^^)), A(^) ^ (4^)) = Tr(p ■ 1^ ® 42)) 
A^;^-* = {^J^-* ® ^-"^ = Tr(p ■ J^^-*), the correlation components. Also by partial 
tracing we define p^^^ = Tr2p, p^^-* = Trip. To proceed with the definition of the 
correlation index we view the space of matrices p G u{Ni) ® u{N2) = ^, as a norm 
space with Hilbert-Schmidt (HS) norm. 



1^1 



(2) 



<A,A> = (TtA^A) 



1/2 



Af2 

, E 



IJ I 



(20) 



for A = (ttij) G Q. This is essentially a Frobenius type matrix norm, which is unitarily 
invariant i.e \\UAY\\ = \\A\\, for U,Y unitary (the lower index of the norm will be 
omitted hereafter). Then we propose the following 

Definition. The correlation scalar index of two coupled finite quantum systems 
been in a mixed state p is defined as 



^^IIApll 



\P-P 



(1) 



P 



{2)||2 



(21) 



Index S provides us with a measure of correlation between the coupled systems in 
terms of the difference of the factorized partial density matrices from the density of 
the composite system. It is cast in the form 



S 



\\P\\ 



,(2)^ 



+ iip«inip(^)i|2 



2Tr(p ■ p^- (X> p 

y y- r;,(l,2) _ ;^(l);^(2) 



.{1,2) _ .(1) ,(2) ■ 

'^Ni-m,N2-n '^Ni-m'^N2-n 



(22) 



mG2 
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where A'^i, A''2-modulo arithmetic apphes in the respective indices. 
The index £ vanishes for product states and by using the reahty conditions of the Aj's 
i.e X^"^* = A^L^^, 1^ = 1,2, A^^;^)* = A^fi^^^^-n we introduce the matrix Amn ■ = 
^mn - = (jm ® 4^^) " (jm) (jn^)^ ^ud re-express the index in the form 



£ = TrAAt . (23) 

This last expression suggests first, that the index £ is determined by the trace of the 
covariance matrix of local observables J^-* and J!^\ and second that it is invariant under 
general unitary tranformations of the group U{Ni- N2), i.e A —>■ WAU ; A''" U'^A'^U, 
with U e U{Ni ■ N2) C U{Ni) (g) U{N2). The last inclusion describes the fact that the 
invariance unitary group of is in general larger than the local unitary transformations 
in which case the symmetry group factorizes (c./ [Q). 

Extensions to three and more coupled systems is straightforward. For three systems 
e.g the composite density matrix involves terms of the operator basis where the J^'s 
are embedded in all possible ways in the 3-tensor space. Also for the reduced matrices 
there are various possibilities in this case i.e p*'-' = Tr^p and = Ttj^p, with cyclic 
permutations of {i,j,k) = (1,2,3). This gives rise to different correlation indices i.e 



£l2S = \\P - P^'^ ^ P^'^^ ® P^'X. 

£m) = \\p-P^'^®P^''^\\h, 

^2(13) = \\P - P''' ^ P^'^h , 



(2) 



^3(12) = \\p-p^'^®pmh- (24) 

Closing we should mention that the correlation index can be expressed in terms 
of the P function of the involved density matrices, associated with the SU{2) group 
coherent state of dimension A^. This possibility as will be explained elsewhere [^, is 
based on the fact the the su{N) algebra generators used here in the expansion of the 
A^-dim density matrices, can be embedded (by means of the polar decomposition of 
the su{2) algebra), into the enveloping algebra U{su{2)). Examples of the finite case 
together with extensions to infinite dimensional quantum systems will also be reported 
elsewhere. 
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Figure captions. 

Fig. 1. Induction of CRT into the Heisenberg group. 
Fig. 2. Component version of the induction of fig.l. 
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